Hard Probe of Geometry and Fluctuations from RHIC to LHC 
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We study the event-by-event hard probe of geometry and fluctuations in the initial condition of 
heavy ion collisions from RHIC to LHC energies. Results for the dominant harmonic, V2 at high p t 
from jet quenching models with varied path-length and matter-density dependences are compared 
with precise data sets at both RHIC and LHC. We find that too strong a path-length dependence 
(e.g. cubic) is ruled out by data at LHC, while the model with strong near-T c enhancement of jet- 
medium interaction fully describes data. Quantification of azimuthal harmonics v n (n=l, 2,3,4,5, 6) 
of high p t hadrons are presented for LHC 2.76 TeV as well as 5.5 TeV collisions. We also point out 
that the near-T c enhancement model naturally expects a less color-opaque medium at LHC. 



Introduction. — Searching for new forms of matter is a 
fundamental quest. In the strong interaction sector of the 
Standard Model described by the Quantum Chromody- 
namics (QCD), various forms of QCD matter may exist 
in the Nature, e.g. inside the compact stars and in the 
early Universe a few microseconds after the "Big Bang" . 
Studying the "condensed matter physics of QCD" has 
been essential for advancing our understanding of mat- 
ter and of the QCD dynamics as well. A highly nontriv- 
ial prediction based on the salient feature of QCD, the 
asymptotic freedom is that we shall expect a new, 
deconfined, and weakly coupled phase of QCD matter at 
asymptotically high temper ature 0, A deconfinement 
transition at certain temperature T c and a quark-gluon 
plasma (QGP) phase above that temperature [3| are ex- 
pected and have been extensively studied using lattice 
QCD simulations. Experimentally such hot QCD matter 
has been created in heavy ion collisions ( "Little Bang" ) , 
with its many properties measured in the past decade at 
the Relativistic Heavy Ion Collider (RHIC) @] and now 
also at the Large Hadron Collider (LHC) [6J. 

The era of RHIC plus LHC provides unique opportu- 
nities for uncovering the underlying structures of hot de- 
conhned QCD matter and deepening our understanding 
of how QCD operates in a strongly interacting many- 
body setting. The various major findings at RHIC (with 
AuAu collisions up to y/s — 200GeV) coherently hint at 
a strongly coupled QGP (sQGP) [7[ above but close to T c 
which is drastically different from the naively expected 
"asymptotically free matter" (AFM). While a full micro- 
scopic picture of the sQGP is yet to come, there have been 
important progresses suggesting that the QCD plasma 
near T c is an emergent matter dominated by dense and 
light (chromo-)magnetic monopoles @. These thermal 
monopoles are perculiar to the region above but close to 
T c , and their condensation at T c marks the onset of con- 
finement. With LHC colliding PbPb at y/s = 2.76 TeV 
now (and yfs = 5.5 TeV in the future), thus creating an 
even hotter QGP of much higher density, it is tempting to 
ask: if the sQGP at RHIC bears peculiar near-T c nature, 
is that quickly turning off with noticeable changes in the 



QGP properties at LHC which shifts to higher T further 
away from T c , and how much is the QGP at LHC get- 
ting closer to the AFM? This paper attempts to extract 
insights into these questions by studying the hard probe 
of geometry and flucutations at both collision energies. 

Jet quenching and geometric tomography. — Just like 
the X-ray imaging of normal materials, the highly en- 
ergetic partons produced in the initial binary hard colli- 
sions provide natural imaging tools to study the created 
hot matter in heavy ion collisions. Such a partonic jet 
carries an energy much higher than the medium energy 
scale and when it travels through the dense medium it 
experiences multiple collisions with medium constituents 
and loses its energy significantly (i.e. jet quenching) [§]. 
As a result the hadrons eventually produced from the jet 
will differ from the case without medium effects e.g. in 
a proton-proton collision at the same energy. Measuring 
such differences is a most useful way of learning about the 
properties of medium and the jet-medium interactions. A 
conventional observable to quantify the jet quenching is 
the nuclear modification factor Raa defined as: 



R A A(pt,(f>,v) 



d 2 N AA /dp t d<j)dri 
= T AA ■ d 2 a NN /dp t d(j)dri 



(1) 



where in the denominator the the nuclear overlap func- 
tion Taa scales up single Nucleon-Nucleon(NN) cross sec- 
tion to Nucleus-Nucleus(AA) according to the expected 
number of binary NN collisions assuming no modifica- 
tion. Thus a value of Raa smaller (larger) than unity 
means suppression (enhancement) due to medium effect. 
For physics related to jet quenching, we focus on Raa 
measured for detected hadrons with large enough trans- 
verse momenta (e.g. pt > 6GeV for RHIC andp t > 8GeV 
for LHC) and for r/ at mid-rapidity (as per most detec- 
tors). The 4> is the azimuthal angle of the measured 
hadron's transverse momentum, and one may examine 
the Raa either integrated or differential in </>. 

Significant suppression of hadron production at large 
transverse momentum p t , was first observed at RHIC, 
with Raa reaching about 0.18 in the most central col- 
lisions. Measurements of charged particle, identified 
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hadron, heavy flavor production, photon production, 
triggered di-hadron correlations have coherently pointed 
to a created medium that is extremely opaque to "col- 
ored" hard probe. For reviews see e.g. [lOj. New and 
extensive LHC data on various hard probe observables 
have shown similar strong jet quenching [llj, and the apt 
question is to quantitatively analyze whether the hotter 
medium becomes less opaque or not (in the context that 
medium density roughly doubles). 

A very powerful idea in jet quenching study is the 

The hot 
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so-called geometric tomography [1^ 
medium created in a heavy ion collision event is gener- 
ally anisotropic in the transverse plane (perpendicular to 
the collision beam axis), therefore high energy partons 
traversing the medium along different azimuthal direc- 
tions will "see" different medium thickness and thus lose 
different amount of energy. That will lead to a measur- 
able anisotropy in the RAA(<f)- It is well-known that the 
dominant geometric anisotropy is the elliptic component 
~ cos 2(0 — ^ 2 ) (with ^2 related to matter anisotropy 
axis) and the coefficient is the v 2 for high p t hadrons. 
A lot of studies |T3 - [2l| have shown that such geometric 
observable of jet quenching is very sensitive to the un- 
derlying dynamics of jet energy loss e.g. its dependence 
on the in-medium path length. More recently the idea 
has been extended to systematically quantify the jet re- 
sponse to various harmonic components of the anisotropy 
arising from the strong event- by- event fluctuations in the 
initial condition of heavy ion collisions [22], [23| . 

Near-T c enhancement of jet quenching. — Substantial 
studies of jet quenching azimuthal anisotropy specified by 
v-i at highpt at RHIC, however, led to a clear discrepancy 
between various model results and the already quite ac- 
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curate PHENIX data extending to ~ 20GeV 
around 2008. All previous jet quenching models, with ei- 
ther linear or the LPM-induced quadratic[25] path-length 
dependence, underproduced the v 2 at highp t significantly 
(often by a factor of 2) after properly constrained in the 
overall opaqueness by (^-integrated Raa- Efforts toward 
reconciling Raa and v% at high p t fostered a more radi- 
cal proposal in (li| that is in contrary to the assumption 
(taken for granted in all previous considerations) that the 
energy loss is simply proportional to plasma constituent 
density (e.g. as per entropy density s). Instead, the key 
insight of [151 ] is that the jet-medium interaction has non- 
trivial dependence on matter density and particularly is 
strongly enhanced in the near-T c region. Such near-T c 
enhancement of jet quenching, in analogy with the well- 
known "critical opalescence" , is well motivated by the 
aforementioned emergence of magnetic monopoles in the 
same regime Q. Phenomenologically this model for the 
first time achieved a simultaneous description of Raa 
and V2 at high p t . It shall be emphasized that the effect 
of near-T c enhancement on jet quenching anisotropy is 
robust and generic: similar successes were reported by 
incorporating the near-T c enhancement in a variety of 



different approaches for jet quenching [261 ]. 

For a more concrete discussion on the path-length de- 
pendence and matter-density dependence of jet quench- 
ing, let us adopt here the geometric models that have 
been widely used and successful in describing the gross 
features of jet quenching [l4T - [l9j . In such a model one 
assumes that the final energy Ef of a jet with initial en- 
ergy Ei after traveling an in-medium path P (specified 
by the jet initial spot and momentum direction) can be 
parameterized as E f — Ei x fp with the fp given by 



fp = exp 



k[s(1)} s{1) l m dl 



(2) 



Here s(l) is the local entropy density along the jet path, 
while the n(s) represents the local jet-medium interac- 
tion strength which as a property of underlying matter 
depends on the local density s(l). We choose to ex- 
plicitly separate out the density s itself, and the com- 
bination k(s) s approximately corresponds to q in many 
jet quenching models. Different choices of m and k(s) 
mean different path-length and matter-density depen- 
dences, and here we consider three classes of models. The 
near-r c enhancement (NTcE) model as in [l5| [l9j is im- 
plemented by assuming m = 1 (i.e. quadratic) and intro- 
ducing a strong jet quenching component in the vicinity 
of T c (with density s c and span of s w ) via 

k(s) = /to[l + £ exp(-(s - s c ) 2 /s 2 w )] , (3) 

with £ = 6, s c = 7//m 3 , and s w = 2/fm 3 . (see [H, H] 
for the details.) For contrast, we also consider two other 
classes of models that both assume n(s) — kq being a 
constant independent of density, while have m = 1 and 
m = 2 respectively, referred to as L 2 and L 3 models here- 
after. The parameter kq controls the overall opaqueness 
in each of the three models and will be fixed by Raa ~ 
0.18 in the 0-5% collisions at RHIC y/s = 200GeV. The 
L 2 model represents the generic feature of most radiative 
energy loss models with LPM effect, while the L 3 model is 
motivated by certain energy loss calculations for strongly 
coupled Yang-Mills plasma based on the AdS/CFT corre- 
spondence [161 127 [ . It is worth mentioning, though, that 
at RHIC it turns out the L 3 model is also able to describe 
v 2 at high p t [l|| HH ■ The reason may be that the L 3 de- 
pendence effectively enhances the later time quenching 
which mimics the similar effect from near-T c enhance- 
ment. Therefore further discrimination between the two 
models is called for and LHC test is crucial. 

Hard probe of geometry and fluctuations from RHIC to 
LHC. — In this study we report the first event- by- event 
quantification of the azimuthal anisotropy in jet quench- 
ing due to both geometry and fluctuations at both RHIC 
and LHC energies. In particular we will show results for 
V2 at high p t from all three geometric models, the NTcE, 
L 2 and L 3 models, and compare with experimental data. 

In general for a given event, there are strong initial 
fluctuations in both the participant density profile (which 
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FIG. 1: (color online) The V2 of high p t hadrons at RHIC 
0.2 TeV (upper) and LHC 2.76 TeV (lower) collisions com- 
puted from NTcE (thick solid blue), L 2 (dashed red) as well as 
L 3 (thin solid black) models. The results are compared with 
various data: PHENIX 6 — 9 (magenta boxes) and > 9 GeV 
(green boxes) for the RHIC (upper); ALICE 10-20 GeV (pur- 
ple asterisks), ATLAS 8-12 GeV and 9-20 GeV (open and 
filled magenta diamonds), as well as CMS 9 and 11 GeV bins 
(open and filled green boxes) for the LHC (lower). 

dominates the soft matter density distribution) and the 
collision density profile (which dominates the jet spot 
distribution). The jet quenching anisotropy comes from 
the convolution of both, and the resulting Raa{4>) can 
be evaluated on an event-by-event basis as follows: 



R AA {cj>) =< (f P J 



(4) 



where < >p , means averaging over all jet paths with 
azimuthal orientation </> and including all the possible 
initial jet production spots (distributed according to the 
binary collision density in the same event). The exponent 
n comes from reference p-p spectrum at the same collision 
energy: n s» 8.1, 6.0, and 5.4 for y/s = 0.2, 2.76, and 5.5 
TeV (see e.g. Estill )■ The so-obtained Raa{4>) in 
each event can be further Fourier decomposed as: 



R AA ((f>) — Raa ( 1 



71=1 



The overall quenching Raa as well as the azimuthal har- 
monics v n and the corresponding n-axis tp^ can then be 
determined from the above in each event, followed by av- 
erage over events ( ~ 10 4 for each impact parameter in 
this study). The second harmonic v 2 is the most robust 
and reflects the hard probe of anisotropy both from ge- 
ometry and fluctuation, while all the other harmonics are 
jet responses to fluctuations and provide further insights 
into the initial conditions in addition to what have been 
learned from the bulk collective expansion dynamics [3lj | . 

To perform the event-by-event simulation, we have 
used the standard Monte-Carlo Glauber model to gener- 
ate fluctuating initial conditions 32] and followed most 
hydrodynamics literature for setting the relevant imple- 
mentation procedures and parameters 33|, |34[ . The NN 
cross section Oi ne x is set as ai ne i = 42, 62, and 66 nib for 
y/s = 0.2, 2.76, and 5.5 TeV. With calculated partici- 
pant density p p (r- L ) and binary collision density p c (r^-) 
in the transverse plane, the initial entropy density at 
equilibrium time ro = 0.6 fm/c (assumed the same for 
different energies) is scaled with p p (r _L ) at RHIC, and 
(1 - S)/2 x pp(r^) + 5 x PcO^) at LHC (S = 0.118 for 
both 2.76 and 5.5 TeV cases) with the proportionality 
constants determined from multiplicities at different en- 
ergies 33]. Energy loss in the pre-equilibrium stage is 
quite possible but quantitatively uncertain, and may be 
improved in the future with developing understanding 
of the thermalization process [13]. We've adopted the 
strategy in [lH ] to "turn on" the pre-equilibrium density 
for t < To with a linear increase and then for t > tq 
to decreases the entropy density as 1 /r thus taking into 
account the boost-invariant longitudinal expansion. 

Let us now first show in Fig. [1] the results for v 2 at 
high p t from all three models and compare them with 

It shall 
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data from both RHIC |28( and LHC 
be emphasized again that all models have their parame- 
ters fixed only with RHIC — 5% data once and for all. 
For the RHIC comparison, the simple L 2 model clearly 
fails to describe data while the NTcE and L 3 models 
are indeed both compatible with data. Moving to the 
LHC comparison, however, one immediately sees that too 
strong a path-length dependence (i.e. the L 3 model) is 
ruled out, while both the NTcE and L 2 models are now 
consistent with data. With RHIC plus LHC, therefore, 
only the model with LPM-type quadratic path-length de- 
pendence and strong near-T c enhancement of jet-medium 
interaction fully describes data at both collision energies. 

We further come to fully quantify other azimuthal har- 
monics w I j(n=l,2,3,4,5,6) of high p t hadrons based on 
our event-by-event computation in NTcE model. The 
results for LHC 2.76 TeV collisions are shown in Fig. [2] 
together with available data for ^3,4,5 HI, 37 1. One shall 
however be cautious with the comparison of higher har- 
monics with data. Different from v 2 for which both the 
final bulk matter event plane and the final hard response 
(quenching) plane are tightly correlated with the initial 
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FIG. 2: (color online) The v n (n=l, 2,3,4,5, 6) of high p t hadrons at LHC 2.76 TeV collisions computed from NTcE model. 
Available data for ^3,4,5 are shown as asterisks (ALICE 10 — 20 GeV), open and filled diamonds (ATLAS 8—12 and 12—16 GeV). 
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FIG. 3: (color online) Predictions from NTcE model for v n 
(n=l,2,3,4,5,6) of high p t hadrons at LHC 5.5 TeV collisions. 

participant plane, for those higher harmonics from fluctu- 
ations the responses of both the soft and hard sectors are 
much less correlated with the initial e n plane and there- 
fore the comparison between computations and event- 
plane-based measurements for higher harmonics could be 
rather tricky [23j, |33[ . Finally in Fig. [3] we show the first 
predictions for v„(n=l,2,3,4,5,6) of high p t hadrons at 
LHC 5.5 TeV collisions, which will be tested in few years. 

Summary and Discussions. — In summary we've stud- 
ied the event- by- event hard probe of the geometry and 
fluctuations in the initial condition of heavy ion collisions. 
With precise data sets for high p t azimuthal anisotropy at 



RHIC plus LHC, jet quenching models with varied path- 
length dependence and matter-density dependence can 
be discriminated. We've found that too strong a path- 
length dependence (e.g. cubic) is ruled out by data at 
LHC, while the model with strong near-T c enhancement 
of jet-medium interaction fully describes data at both 
collision energies. Finally full quantification of the az- 
imuthal harmonics v n (n=l,2,3,4,5,6) of high p t hadrons 
are presented for LHC 2.76 TeV as well as 5.5 TeV col- 
lisions. Admittedly the modeling in this study is not 
the most sophisticated, but we expect the gross conclu- 
sion from studying generic geometric features to be ro- 
bust. Further improvements (e.g. including energy de- 
pendence, integration with realistic hydrodynamic mod- 
eling, possibility of incorporating the NTcE component 
with other widely used jet quenching schemes, etc) are 
underway and will be reported in the future. 

The near-T c enhancement model implies a strong de- 
crease of jet-medium interaction with hotter tempera- 
ture and naturally expects a less color-opaque medium at 
LHC despite only modest increase of temperature from 
RHIC UM- Consistent messages have been reported re- 
cently from a variety of independent jet quenching studies 
[39M42| . The underlying picture of emergent plasma near 
T c may imply a more rapid running than what would be 
expected from perturbative picture and it might be pos- 
sible that in LHC top energy heavy ion runs the created 
quark-gluon plasma could be considerably closer toward 
the long expected asymptotically free matter. 
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